Controlled wave-packet manipulation with driven optical lattices 
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Motivated by recent experimental progress achieved witii ultracold atoms in kiloiiertz-driven op- 
tical lattices, we provide a theoretical discussion of mechanisms governing the response of a particle 
in a cosine lattice potential to strong forcing pulses with smooth envelope. Such pulses effectuate 
adiabatic motion of a wave packet's momentum distribution on quasienergy surfaces created by 
spatiotemporal Bloch waves. Deviations from adiabaticity can then deliberately be exploited for 
exerting coherent control and for reaching target states which may not be accessible by other means. 
As one particular example, we consider an analog of the vr pulses known from optical resonance. We 
also suggest adapting further techniques previously developed for controlling atomic and molecular 
dynamics by laser pulses to the coherent control of matter waves in shaken optical lattices. 

PACS numbers: 67.85. Hj, 42. 50. Hz, 32.80.Qk, 03.75.Lm 
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I. INTRODUCTION 

Boosted by the seminal observation of the quantum 
phase transition from a superfluid to a Mott insulator in 
a gas of ultracold ^''Rb atoms trapped by an optical lat- 
tice potential [ij , the experimental and theoretical study 
of ultracold atoms in optical lattices has matured into a 
major area of contemporary research 043 ■ To a large 
extent, this field is driven by the promise of simulating 
complex condensed-matter systems and obtaining novel 
insight into phenomena which hitherto are not under- 
stood, such as high-temperature superconductivity. 

At present, evidence is accumulating which suggests 
that this field is developing a new branch, aiming at the 
coherent control of mesoscopic matter waves in optical 
lattices through the application of time-periodic forces, 
with driving frequencies in the lower kilohertz regime. 
While it had been pointed out earlier that a metal- 
insulator-like transition undergone by ultracold atoms 
in quasiperiodic optical lattices should be controllable by 
adjusting the amplitude of a sinusoidal drive {5|, experi- 
mental work in this direction increased pace only in 2007, 
with the clear-cut observation of forcing- induced dynami- 
cal suppression of tunneling, and even reversal of the sign 
of the tunnelingmatrix element, by Arimondo, Morsch, 
and co-workers ^0|- This group also has documented an 
analog of photon-assisted tunneling with Bose-Einstein 
condensates in shaken optical lattices Q, and has veri- 
fied that the superfluid-to-Mott-insulator transition can 
be coherently controlled by suitably "dressing" a matter 
wave in an optical lattice B , taking up a theoretical pro- 
posal by Eckardt et al. [10]. The very same principle un- 
derlying this form of coherent control has quite recently 
been exploited successfully for emulating frustrated mag- 
netism in driven triangular optical lattices j7l] |. which 
may well be regarded as a guiding landmark example 
of quantum simulation. Moreover, there now exist first 
experimental results demonstrating active control of cor- 
related tunneling in ac-driven optical lattices [l2;] • 

These experimental advances concerning time- 
periodically-driven matter waves are accompanied by 



growing theoretical efforts. For example. Kudo et al. 
have investigated the possibility of driving-induced 
control of bound-pair transport [Tsj , while Tokuno 
and Giamarchi have studied a kind of spectroscopy for 
cold atoms in periodically phase-modulated optical lat- 
tices [lil- Moreover, Tsuji et al. have pointed out that 
ac forcing may even change the interparticle interaction 
from repulsive to attractive, possibly allowing one to 
simulate an effectively attractive Hubbard model with 
a temperature below the superconducting transition 
temperature So far, all these considerations merely 
involve strict ac forcing with a constant amplitude. In 
analogy to the physics of atoms and molecules interact- 
ing with laser pulses, here we suggest that many more 
control options should become available when ultracold 
atoms in optical lattices are subjected to forcing pulses 
with a deliberately shaped envelope. 

Such attempts to gain coherent control over meso- 
scopic matter waves call for a systematic theory of the 
response of ultracold many-body systems to nonpertur- 
batively strong external forcing. Although one may ob- 
tain some insight from drastically simplified model sys- 
tems [l^, and substantial progress is being made now 
with the help of advanced numerical schemes this 
goal still is far from accomplishment. In this situation, 
an intermediate step suggests itself: In experiments with 
sufficiently dilute Bose-Einstein condensates in optical 
lattices, or with condensates for which the interparticle 
s-wave scattering length has been tuned close to zero by 
means of a Feshbach resonance, one may ignore inter- 
action effects altogether, and can observe typical single- 
particle phenomena, such as ordinary or "super" Bloch 
oscillations, with condensates jl8n21j . Thus, it seems ad- 
visable to undertake a comprehensive theoretical study 
of the possibilities of coherent control of single-particle 
dynamics in forced optical lattices. The results of such 
a study can then immediately be tested in "interaction- 
free" condensate experiments and may help one to dis- 
entangle genuine many-body effects at a later stage. 

This is the step we are going to take in the present 
paper. Building on our previous work p^ . here we pro- 
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vide a detailed picture of basic mechanisms which im- 
ply single-particle state control: We first demonstrate in 
Sec. ini the feasibility of adiabatic transport of momen- 
tum distributions on quasienergy surfaces corresponding 
to time-periodically-forced optical lattices. This option 
is opened up by the existence of a basis of spatiotempo- 
ral Bloch waves, that is, of Bloch-like states which em- 
body both the spatial periodicity of the lattice and the 
temporal periodicity of a driving force on equal footing; 
such states constitute the foundation of our analysis [23 . 
We then establish in Sec. IIIII an analog of the tt pulses 
known from the theory of optical resonance 2J] and out- 
line how to utilize avoided crossings of quasienergy sur- 
faces for "cutting out" parts of an initially given mo- 
mentum distribution. Section ITVl brief! v addresses effects 
connected to the phase of the driving force. Taken to- 
gether, our findings indicate that there is a high potential 
for transferring well-established methods currently used 
for manipulating and controlling atoms and molecules 
by specifically designed laser pulses (25l - [29j to the newly 
emerging field of manipulating and coherently controlling 
mesoscopic matter waves in optical lattices by specifically 
tailored forces; this prospect is put forward in our con- 
clusions. 



II. ADIABATIC TRANSPORT OF MOMENTUM 
DISTRIBUTIONS 

The starting point of our considerations is a single par- 
ticle of mass m moving in a one-dimensional optical lat- 
tice potential 0-11] 
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cos {2ki^x) , 



(1) 



where the lattice depth Vq is proportional to the inten- 
sity of the laser radiation generating the lattice, and /cl 
denotes the corresponding wave number. Thus, the po- 
tential is periodic with lattice constant a = n/k-i^, so that 
V{x) = V{x + a). Moreover, the particle is subjected to 
a spatially homogeneous inertial force F{t), which can 
be applied by accelerating the lattice in the laboratory 
frame [s^ . After transforming to a frame of reference co- 
moving with the lattice, the Hamiltonian of the system 
is given by 



H{x,t) = ^ + V{x) - F{t)x 
2m 



(2) 



If we denote the solution to the Schrodinger equation 
pertaining to Eq. ^ by •il}{x,t) and perform the unitary 
transformation 



tpix, t) — exp 



ATF{T)]i>{x,t) , (3) 



the transformed functions iplxjt) obey a Schrodinger 
equation with the new Hamiltonian 

H{x,t) = ^(p+ fdrFir] 



) +Vix). (4) 



The traditional solid-state approach to monitoring the 
wave-packet dynamics now is as follows: The unforced 
lattice possesses improper energy eigenstates Xn,kix) 
which have the form of Bloch waves [3TI - [33| . that is, 
of plane waves which are modulated by lattice-periodic 
functions w„,/c(x) — w„,/c(x -I- a), so that 



Akx 



Vn.k{x) 



(5) 



these waves solve the time-independent Schrodinger 
equation 



h V{x) ] Xn.k{x) = En{k)Xn.k{x) ■ (6) 

2m 



Here n is a band index and k a wave number, so that 
En{k) is the energy dispersion relation of the nth Bloch 
band. Owing to the periodicity of the lattice, the wave 
numbers can be restricted to the first quasimomentum 
Brillouin zone B ~ [— vr/a, +7r/a[. In addition, we require 
the normalization 



'^^XXn'k'{x)Xn,k{x) = 5n,n'5{k - k') . (7) 



When an arbitrary given wave packet tp{x,t) is expanded 
with respect to the Bloch basis in the form 



i'{x,t) = j^'^^9n{k,t)Xn,k{^ 



(8) 



this convention ([7]) makes sure that the momentum dis- 
tributions |(7^(fc,i)p are normalized according to 



I dfc|g«(fc,t)p = l. 



(9) 



In solid-state physics, the expansion (|8]) is known as 
the crystal-momentum representation of the wave packet 
^j{x,t). 

When there is no external forcing, F{t) = 0, the 
time dependence of the expansion coefficients g^{k,t) in 
Eq. ([U simply reads 



(10) 



For studying the dynamics under the action of a force 
F(t), let us at this point assume that the wave packet ini- 
tially occupies only one band with a particular index n, 
and that the force remains so weak that it does not in- 
duce substantial interband transitions. Then Bloch's ac- 
celeration theorem [sH - lssj comes into play: The packet's 
center wave number in k space, given by the first moment 



fce(i)= / dkk\g^{k,t)\' 

JB 



(11) 



then evolves according to the semiclassical law 

hk,{t) = F{t) . (12) 
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For example, a constant force leads to a linearly increas- 
ing kc(t), which, in its turn, gives rise to Bloch oscilla- 
tions in real space. 

Although this time-honored approach has many 
virtues, for our purposes it is advantageous to look at the 
wave-packet dynamics from a different angle. In view of 
the goal to exert coherent control on the lattice atom, it 
is quite natural to specifically consider sinusoidal forces 
F{t) = sin(wt) in the first place, since then the Hamil- 
tonian ^ is of the familiar form which also describes 
a charged particle in a monochromatic classical radia- 
tion field within the dipole approximation. The most 
conspicuous difference concerns the frequencies: Typi- 
cal frequencies for driving optical lattices [6l-[T2j fall into 
the lower kilohertz regime, about 11 orders of magni- 
tude lower than optical frequencies. In the following, we 
merely require that the force be periodic in time with 
period T, so that F{t) — F{t + T), and we assume that 
its one-cycle average vanishes, so that 

^j\tF{t)=Q. (13) 

With these specifications, the transformed Hamilto- 
nian (j4]) is periodic in space as well as in time, H{x, t) = 
H{x + a,t) = H{x,t + T). While the ordinary Bloch 
waves (O account for the spatial periodicity only, the 
mathematical Floquet theorem governing the structure 
of solutions to differential equations with periodic co- 
efficients 34 36] can now be invoked to simultaneously 
incorporate both the spatial and the temporal periodic- 
ity, resulting in a set of solutions to the time-dependent 
Schrodinger equation of the suggestive form 

i>n,k{^^ = exp [ifca; - kn{k)t/h] u„,k{x, t) (14) 

with biperiodic functions Un,k{x,t) which reflect the two 
translational symmetries, Un,k{x,t) = u„.fe(x -I- a, t) = 
Un,k{x, t + T). We refer to these solutions as spatio- 
temporal Bloch waves [23j . The quantities £n{k) deter- 
mining the linear growth of the phase factors with time 
are commonly known as quasienergies (stI . [38j . They are 
obtained by solving the eigenvalue problem 

d \ 

H{x,t) - j fn,kix,t) = en{k)^pn,kix,t) , (15) 

where the functions ipn,k{x,t) = exp (ikx) Un,k{x,t) 
denote the spatial parts of the spatiotemporal Bloch 
waves (fT4|) . Fully in accordance with our rationale, this 
eigenvalue problem (jlSp is posed in an extended Hilbert 
space which puts position x and time t on equal foot- 
ing [3^. Because of the periodicity of H{x,t) in time, 
the eigenvalues £n(fc) are defined up to an integer multi- 
ple of huj, with u! — 27r/T, which means that there also 
is a Brillouin-zone scheme for the quasienergies, with the 
fundamental zone Q = [—huj/2, +hLu/2[, in analogy to the 
fundamental quasimomentum zone B = [— 7r/a, +Tr/a[. 
Evidently, Eq. (IT5|) now takes the place of the traditional 



eigenvalue equation ([6]), and the spatiotemporal Bloch 
waves replace the ordinary Bloch waves ([5]). Con- 
sequently, we abandon the standard crystal-momentum 
representation and instead perform expansions of 
given wave packets ijj{x,t) in this new basis: Fixing, in 
analogy to the previous Eq. ([7]), the normalization 




dxifl, ,^,{x,t)tpn,k{x,t) = Sn,n'S {k - k') , (16) 

' n 



we thus arrive at the Floquet representation [23| 

'^ix,t) = ^ J dkgn{k,t)(p„^k{x,t) . (17) 

When the amplitude of the driving force goes to zero, 
the functions ifin,k{x,t) reduce to the Bloch waves ([5|), 
and the quasienergies Snik) approach the energies i?„(fc), 
modulo huj. Therefore, in this limit the Bloch expan- 
sion dHl) coincides with the Floquet expansion p7)) . How- 
ever, in the presence of a strictly T-periodic force, such 
as F(t) = i^ac sin(wt), we now have two different pictures 
of the same wave-packet dynamics: Within the crystal- 
momentum approach, a single-band wave packet is de- 
scribed by a momentum distribution \g^{k,t)\'^; the cen- 
ter of this distribution moves in k space according to the 
acceleration theorem (IT^ . In contrast, within the Flo- 
quet picture one merely has 

g„(fc,i)=gn(fc,0)e--"W*/\ (18) 

so that the Floquet distribution \gn{k,t)\'^ = |f;„(fc,0)p 
does not move at all, but stays perfectly constant in time, 
with the response to the oscillating force already being 
incorporated into the basis states ((T4)) . 

Clearly, both approaches are mathematically equiva- 
lent. But it is the second one which allows us to make 
further contact with advanced techniques developed for 
studying the interaction of atoms and molecules with 
laser radiation, and for developing schemes for coherent 
control of ultracold atoms in driven optical lattices. 

Such schemes naturally will involve pulses of driving 
forces, that is, non-periodic forcing. As a simple case, we 
may consider pulses of the form 

F(0=Far^s(t)sin(c.i), (19) 

where the dimensionless shape function s{t) vanishes be- 
fore and after the pulse, s{t) — for both t < and 
t > Tp, say, and is normalized such that its maximum 
value is 1, implying that F^^^ is the maximum amplitude 
encountered during the pulse. 

It is then of key importance to note that the Floquet 
picture is meaningful not only for perfectly time-periodic 
forces F{t) = F{t + T), but also for situations in which 
one or more system parameters change slowly, that is, un- 
dergo only minor variations during one cycle T [40j . This 
is the case, for instance, if the pulse F{t) is equipped with 
a "slowly" varying envelope s(t). Then one considers not 
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only one single eigenvalue problem (jlSp corresponding to 
one particular amplitude i^acj but rather the family of all 
such eigenvalue problems with < Fac < FJ^^^. This 
gives a basis of spatiotemporal Bloch waves for each in- 
stantaneous value of the amplitude; taken together, these 
bases serve as a "moving frame of reference" with respect 
to which the wave packet can evolve adiabatically: Under 
pulse conditions enabling adiabaticity, the Floquet mo- 
mentum distributions \gnik,t)\'^ remain almost constant 
in time, provided the expansion (|17p refers at each mo- 
ment to that basis of spatiotemporal Bloch waves which 
is obtained by fixing the slowly varying amplitude at 
its momentary value. Pictorially speaking, the instanta- 
neous eigenvalues e^"" (fc), considered as functions of wave 
number k and driving amplitude i^ac, form quasienergy 
surfaces on which the momentum distribution can move 
almost without change of shape in response to smooth 
variations of the envelope s(t). 

To see what this means in practice, we consider a co- 
sine lattice ([T|) with depth Vq = 5.7 E,-, where R- = 
h^k'^/{2m) denotes the single-photon recoil energy [2|tJ]; 
this depth is routinely being realized in current experi- 
ments 0-[i[ni,[il- The width of the lowest Bloch band 
then amounts to Ei{ki,) - Ei{Q) = 0.220 E'r, while the 
lowest band gap is i?2 (fcL) — (^l) = 2.763 E^. The max- 
imum separation of the lowest two bands, encountered 
in the Brillouin zone center, figures as -£2(0) — i?i(0) — 
4.690 i?r [12] ■ We then take a wave packet prepared at 
time i = in the lowest Bloch band n = 1, 

^{x,0) = ^ J^dkgfik,0)xiAx) , (20) 
with a Gaussian initial momentum distribution 

gf{k,0) = (V^Afc)'^/'exp (21) 

centered around fcc(0)/fcL ~ with width Afc/fcL = 0.1. 
This wave packet is subjected to pulses (ITOl) of length Tp 
with a squared-sine envelope, 

s{t) = sin^ {^T'^ ' ^ - * - ■ (22) 

The frequency selected for the model calculation dis- 
cussed in the following is w = 1.640 well below 
the lowest band gap, while the maximum scaled driving 
amplitude 

pmax_ 

is set to K™'^'^ = 0.8. The pulse length is fixed at 
Tp = 10 T, so that the peak driving strength is reached 
within no more than five cycles. We then compute, on 
the one hand, the momentum distributions jgf (fc,t)p in 
the crystal-momentum representation, and also perform 
the Floquet expansions 

^) = E \l^j/^3n{k, ^'^^nf^ (24) 




FIG. 1: (Color online) Response of the initial wave packet (|20p 
with momentum distribution (|21|l in an optical lattice with 
depth Vo/Er = 5.7 to a short pulse (|19p with nonresonant 
scaled frequency hu/Er — 1.640, maximum scaled amplitude 
^max ^ Q p^jgg length Tp/T = 10. (a) shows the den- 

sity \gf(k,t)\'^ in the crystal-momentum representation. For 
comparison, the white-dashed line is the first moment kc{t), 
as predicted by the acceleration theorem (|12p . (b) depicts the 
Floquet density \gi{k,t)\'^, obtained by expanding the same 
wave packet with respect to the instantaneous spatiotemporal 
Bloch waves. 

with respect to the instantaneous solutions (^^"^ (x, t) to 
the quasienergy equation (fT5)) in order to obtain the cor- 
responding Floquet distributions \gi{k, t)\'^, on the other. 

Figure [1] juxtaposes the results of the two approaches. 
In Fig. [Tfa) we show the evolution of the crystal- 
momentum density |gf (fc,t)p under the pulse. Because 
interband transitions remain negligible for the param- 
eters chosen, the packet's center moves in perfect ac- 
cordance with Bloch's acceleration theorem (|12p . On 
the other hand. Fig. [ijb) depicts the evolution of the 
Floquet density \gi{k,t)\'^. This density remains practi- 
cally constant, indicating almost perfect adiabatic follow- 
ing of ip^Xjt) with respect to the spatiotemporal Bloch 
waves: Despite the short duration of the pulse, the initial 
distribution merely makes an adiabatic excursion on its 
quasienergy surface, returning more or less unaltered. 

Still, Fig.[T]is no more than a look at the same dynam- 
ics from two different viewpoints, and so far neither of 
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Zi^max 



FIG. 2: (Color online) Final escape probabilities from the 
lowest Bloch band of an optical lattice with depth Vo/Er = 
5.7, calculated for the same initial wave packet as considered 
in Fig. [T] after pulses with squared-sine envelope (|22p and 
length Tp/T — 50. Observe the window of almost adiabatic 
response appearing between the two-photon and the single- 
photon resonances. 



these is better than the other. But now comes the cru- 
cial step: Control is exerted by utilizing interband tran- 
sitions. While such transitions fall outside the scope of 
the semiclassical acceleration theorem, which explicitly 
requires a single-band setting, they can be monitored as 
deviations from adiabaticity, caused by near-degeneracies 
of quasienergy surfaces, within the Floquet approach. In 
order to locate the parameters for which such deviations 
occur. Fig. [5] shows the final escape probability from the 
lowest Bloch band in the amplitude- frequency plane [23 | , 
as resulting from the same initial wave packet as con- 
structed above after pulses with the envelope ([2^ . with 
greater length Tp = 50 T. Most notably, the "single- 
photon resonance" with huj = -©2(0) — Ei{0) shows up 
already for quite small driving amplitudes around huj « 
4.690 ii^r; and one observes a sequence of multiphoton- 
like resonances at lower frequencies. Interestingly, there 
also is a pronounced frequency window between the two- 
photon resonance and the single-photon peak which al- 
lows for adiabatic response even to fairly strong pulses 
with K™^^ > 3. This window appears to be most suit- 
able for studying single-band phenomena associated with 
strong forcing, such as the driving- induced reversal of the 
sign of the effective hopping matrix element 6, 7]. 

The Floquet approach now enables one to look into the 
transition dynamics in great detail, and thus to under- 
stand basic principles allowing one to deliberately manip- 
ulate the momentum distribution and to create certain 
desired target states; this will be elaborated in the follow- 
ing section. Before closing the present section, we would 
like to draw an interesting comparison: Aside from de- 
generacies, the adiabatic motion of a wave packet's mo- 
mentum distribution on its quasienergy surface, as vi- 
sualized in Fig. [IJb), seems to resemble the adiabatic 
evolution of molecular states on their Born-Oppenheimer 



potential energy surfaces 4l|. There is, however, an im- 
portant difference: In the case of cold atoms in driven 
optical lattices the concept of adiabatic following [H, HI] 
has to be applied to each wave number k in parallel, 
each one labeling a different spatiotemporal Bloch wave. 
Thus, here we are confronted not with adiabatic follow- 
ing of individual states, but rather with that of a density 
associated with a continuum of quasienergy eigenstates. 



III. TAILORING THE MOMENTUM 
DISTRIBUTION 

Non-adiabatic transitions, which prevent an initial 
momentum distribution from returning practically un- 
changed after a pulse, result from near-degeneracies of 
quasienergy surfaces. Two different cases have to be 
distinguished: Either the near-degeneracy is induced al- 
ready at small driving amplitudes by selecting a resonant 
frequency, or it shows up only under strong nonresonant 
driving, when the ac Stark shift forces two surfaces into 
an avoided crossing (2^. In this section we show that 
either of these scenarios can be exploited for controlling 
and reshaping the k space distribution coherently. In 
all model calculations we consider an optical lattice with 
depth Vq = 5.7 Er and start from the initial Bloch wave 
packet ((20|) with Gaussian coefhcients (|211) , again setting 
Afc/fcL = 0.1. 



A. Resonant forcing 

We now adjust the driving frequency such that hu; — 
£2(0) — Ei{0), so that the lowest two bands are coupled 
resonantly in the center of the quasimomentum Brillouin 
zone, at fc/fcL = 0. The length of the pulses with squared- 
sine envelope is Tp — 50 T. Figure |3] shows the 
resulting final distributions \gf{k,Tp)\'^ and |g^(fc,Tp)p 
for the lowest and for the first excited Bloch band, re- 
spectively, in dependence on the maximum driving am- 
plitude K™^^. One observes a smooth, oscillating exci- 
tation pattern, the first indication of which was already 
visible in Fig. [21 In particular, for K^^'^'^ = 0.186, which 
is the lowest peak amplitude leading to maximum popu- 
lation transfer to the band n = 2, the final excited-band 
distribution is substantially narrower than the original 
one, while the distribution remaining in the lowest band 
is bimodal, corresponding to a wave packet moving in 
two opposite directions. 

With the help of the tools assembled in the preced- 
ing section, these results can be understood in an almost 
intuitive manner, without the need to invoke much for- 
malism. Figure lU^a) depicts the two quasienergy surfaces 
involved in the dynamics, emerging from the two unper- 
turbed Bloch bands n = 1 and n = 2. Because of the 
Brillouin-zone structure of the quasienergy spectrum, the 
"one-photon resonant" driving frequency causes a degen- 
eracy of both surfaces for vanishing instantaneous scaled 
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FIG. 3: (Color online) Final momentum distributions in the 
lowest (a) and in the first excited (b) Bloch band after the 
initial state has been exposed to pulses with squared-sine en- 
velope (I22p and length Tp/T = 50, with varying maximum 
scaled amplitudes < K^'^^ < 0.8. Here the scaled driving 
frequency is hw/Er — 4.690, implying Hw — i?2(0) — £'i(0), 
so that both bands are exactly resonant at k/ki, = 0. The 
lowest panel (c) compares the final distributions in the first 
(dotted) and in the second (dashed) bands to the initial distri- 
bution (full line), for isT™'"' = 0.186. This particular situation 
corresponds to a "tt pulse". 



amplitude Kg = F^cU / (huj) . As a consequence, the initial 
wave packet is not placed on an individual quasienergy 
surface under the action of a pulse, but rather placed 
coherently on both surfaces: A Floquet expansion PT)) 
yields contributions for both n = 1 and n = 2. Both 
parts of the wave function then react adiabatically to 
the slowly varying amplitude, each one picking up its 
own dynamical phase factors, given by the time integrals 
over the instantaneous quasienergies e^'^" {k) encountered. 
When the driving amplitude goes to zero at the end of 
the pulse, both parts of the wave function produce an in- 
terference pattern which determines the final excitation 
probability: For each wave number k sufficiently close to 
resonance, the transition probability to the first excited 




FIG. 4: (Color online) (a) Quasienergy surfaces underlying 
the excitation pattern observed in Fig. [S] The upper sur- 
face originates from the unperturbed Bloch band n = 1, the 
lower one from the band n = 2. Because of the resonant 
frequency, both surfaces are degenerate at fc/ZcL = for van- 
ishing instantaneous amplitude Kq. The quasienergy lines at 
k/kh = are emphasized for better visibility, (b) Section 
through the surfaces at k/ki, — 0, showing the removal of the 
initial degeneracy. 



band is proportional to the expression [44 



1^2 




(k) 



(fc) 



(25) 

As seen in Fig.|Hb), for k/k]^ — the quasienergy differ- 



ence e^ 



10) 



(0) increases linearly with the driving 



amplitude, as is typical for a single-photon resonance '441 . 
Maximum excitation then is obtained when the argument 
of the squared sine in Eq. (1^ equals an odd- integer mul- 
tiple of 7r/2, the first such maximum showing up for 



dt 



Fa, 
£1 



it) 



(k) 



-Fa, 
£2 



(*) 



(fc) 



(26) 



This is reminiscent of the familiar 7r-pulse condition; in- 
deed, when the quasienergies are calculated analytically 
within the rotating- wave approximation, Eq. (|26p reduces 
to the customary area theorem known from optical reso- 
nance [13, Ell • But here we are confronted with the fact 
that this condition (j26l) can not be met simultaneously 
for all components k with one single pulse shape: When 
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FIG. 5: (Color online) Floquet representation of the evolu- 
tion of a wave packet initially prepared in the lowest Bloch 
band, under the action of a resonant tt pulse with K^^^ — 
0.186. This figure shows how the final bimodal distribution 
|(;i(fc,Tp)p depicted in Fig. |3Ic) appears after an initial re- 
duction of the original density, corresponding to the partial 
occupation of the other resonantly coupled quasienergy sur- 
face; and after a period of almost adiabatic motion during the 
middle of the pulse. 



it is satisfied for k/kj^ — 0, the other components of the 
wave packet experience slightly or even strongly different 
quasienergies, depending on its initial width in k space, 
as becomes evident when looking at Fig. Hl^a) . This is ex- 
actly what allows one to "cut out" a part of the momen- 
tum distribution, as was demonstrated in Fig. EJc) : Here 
the pulse shape is such that Eq. (|26l) indeed is satisfied 
for fc/fcL = 0, leading to maximum transition probability 
in the center of the Brillouin zone. In contrast, the ini- 
tial degeneracy at k/ki^ — has no effect on the wings 
of the initial distribution, so that these wings return adi- 
abatically. As a result, the pulse transfers a relatively 
narrow central part of the initial distribution to the first 
excited Bloch band, leaving behind a symmetric bimodal 
distribution in the lowest one. 

It is then of particular interest to monitor the dynamics 
during such a pulse in the Floquet representation, instead 
of merely looking at the final distributions in the usual 
crystal-momentum representation. One such example, 
visualizing the the action of the very "tt pulse" consid- 
ered above, is shown in Fig.[5l Observe how the evolution 
of |(?i(fc, <)P embodies the elements discussed before: The 
distribution soon is reduced to half its initial height, re- 
flecting the occupation of the other quasienergy surface 
at the beginning of the pulse; then stays about constant 
during the pulses' middle part, reflecting approximately 
adiabatic motion; and develops the bimodal pattern only 
at its end, reflecting the final interference. 




^ -0.2 0.2 

FIG. 6: (Color online) Final momentum distributions in the 
lowest (a) and in the first excited (b) Bloch band after the 
initial state has been exposed to pulses with squared-sine en- 
velope (|22p and length Tp/T = 50, with varying maximum 
scaled amplitudes 0.7 < if™^^ < 1.3. Here the nonresonant 
scaled driving frequency is hw/Er — 1.640, as in Fig. [T] (c) 
compares the final distributions in the first (dotted) and in 
the second (dashed) bands to the initial distribution (full line) 
for ifi""" = 1.3. 



B. Nonresonant, strong forcing 

For the following second example of wave packet ma- 
nipulation we again select the driving frequency uj — 
1.640 E^/h, as in the previous calculations having led to 
Fig. [U but now we also consider pulses with larger scaled 
amplitudes K^^^. In Fig. [5] we display the final distri- 
butions \gf{k,Tp)\'^ and |<7^(fc,7p)P as resulting from 
pulses with 0.7 < if^^'' < 1.3. Even for A'J"^'' = 0.8 
the initial distribution returns still undistorted, but for 
^max ^ Q g g^j-Qng iutcrband transitions set in, leading to 
a trimodal excited-band distribution when K^^^ = 1.3. 

Once more the explanation for this response behavior 
is provided by the morphology of the quasienergy sur- 
faces, shown in Fig. [71 The surface with comparatively 
low curvature originates from the Bloch band n = 1 ; this 




FIG. 7: (Color online) (a) Quasienergy surfaces underlying 
the excitation pattern observed in Fig. [S] The almost flat 
surface originates from the comparatively narrow Bloch band 
n = 1, the one with larger curvature from the band n = 2. 
Both surfaces undergo an avoided crossing along a parabola 
with apex ai k/ki, — and Ko ~ 0.9. The quasienergy lines 
at fc/^L = are emphasized for better visibility, (b) Section 
through the surfaces at fc/^L =0, showing the narrow avoided 
crossing. 

surface is penetrated by the one emerging from the Bloch 
band n = 2 along a parabola-shaped line with apex at 
k/ki^ = and Kq « 0.9. Along this line the two surfaces 
exhibit a narrow avoided crossing. The quasienergy rep- 
resentatives plotted in Fig. [7] are shifted against those 
continuously connected to the original energy bands by 
+hu) (n = 1) and by —2huj (n = 2), respectively, so that 
the anticrossing marks a "three-photon resonance." As 
long as the maximum pulse amplitude does not reach the 
apex of the anticrossing parabola, the momentum dis- 
tribution merely moves adiabatically on its quasienergy 
surface, as already demonstrated in Fig. [TJ and returns 
without notable modification. But when K™'^^ > 0.9 
the distribution has to pass the avoided-crossing line, re- 
sulting in partial Landau-Zener transitions to the other 
surface. Thus, the onset of excited-band population in 
the K^^^-k/k]^ plane, as observed in Fig. |6l[b), precisely 
reflects the locus of the band intersection. 

Considering a pulse with K^^^ > 0.9, the different 
k components thus encounter "their" respective avoided 
crossing at different amplitudes, and hence at different 
times. Therefore, the different components acquire quite 
different dynamical phase factors between their first pas- 
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FIG. 8: (Color online) Floquet representation of the evolution 
of a wave packet initially prepared in the lowest Bloch band, 
under the action of a nonresonant pulse with ii'™^" = 1.3. Af- 
ter an initial period of almost adiabatic motion, the center of 
the distribution undergoes partial Landau-Zener transitions 
to the anticrossing surface depicted in Fig.[71 and later returns 
to the initial surface, subject to Stiickelberg oscillations. 

sage through an avoided crossing during the rise of the 
pulse and the second passage taking place during the 
switch-off. Thus, for each k one finds Stiickelberg oscil- 
lations [45|] due to the interference of the parts having 
evolved on the two different surfaces, but their phases 
vary strongly within the Brillouin zone. This feature is 
the origin of the trimodal distribution \gf{k, Tp)p shown 
in Fig.ini^c). As can be clearly seen in Fig.lHl^b), the center 
lobe of this distribution already corresponds to the sec- 
ond Stiickelberg maximum, whereas the two outer lobes 
still are associated with the first. 

In Fig. [8] we show the evolution of the Floquet distribu- 
tion |(7i(fc, <)p during the pulse with maximum amplitude 
^-max — 2.3. Evidently the outer wings of this distribu- 
tion move adiabatically, not encountering the avoided- 
crossing line, whereas the central part jumps to the anti- 
crossing surface with high probability when passing this 
line, returning later when the line is hit a second time. 

We remark that by suitable choices of the frequency 
one may likewise design pulses which involve higher 
quasienergy bands: If the Bloch band ti = 1 is slightly de- 
tuned from a multiphoton resonance with a higher band, 
and if the ac Stark shift forces the quasienergy surfaces 
emerging from these two bands into an anticrossing, one 
can exploit the corresponding multiphotonlike Landau- 
Zener transitions in the same manner as in the example 
considered here. 



IV. PHASE EFFECTS 

The solutions to the quasienergy eigenvalue equa- 
tion (jlSp refer to perfectly periodic driving, and thus do 
not change when the force F{t) = i^ac sm{ujt) is replaced 
by F{t) = i^ac sin(a;t -|- ip), except for a trivial shift of 
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FIG. 9: (Color online) Contour plot of the real-space den- 
sity associated with the lowest Bloch energy band 
(a), and of the density \'\l>i(x ,ty\^ associated with the first ex- 
cited band (b), as resulting from a pulse with hw/E^ — 1.640, 
J^™'"' = 1.5, and T/Tp = 30. The phase ip in Eq. ((27|) has 
been set to zero, (c) shows the gradual loss of population 
from the lowest band. 



the time coordinate. Accordingly, as long as the pulse 
dynamics are fully adiabatic they are not affected by the 
phase if. This is different, however, under non-adiabatic 
conditions. Then the transitions effectuated by a pulse 
may strongly depend on tp, so that this phase offers an 
additional handle of control. For demonstration, we re- 
place the previous pulses ([TO]) by 



F(i) = F--s(i)sinM + vp) 



(27) 



maintaining the interval < i < Tp as the active 
pulse period and employing the same squared-sine en- 
velope (221) as before. Having set oj = 1.640 Sr/ft, 
i^Tjnax ^ 2.5, and T = 30 Tp, Fig. H shows contour plots 
of the densities \ipi{x,t)\'^ and \^p2{x,t)\'^ associated with 
the lowest two Bloch bands, together with the transition 
dynamics, for ip ~ 0. Evidently the transfer of probabil- 
ity to the excited band does not proceed symmetrically 
in space. 

Figure [TO] then depicts the corresponding results for 
(f = TT. While the distribution \4'2{x,t)\'^ here is the 
exact mirror image of that displayed in Fig. |9Ub), the 
depopulation of the lowest band proceeds in exactly the 
same manner as before. This finding is easily explained: 
Replacing ip by p + tt, and simultaneously replacing x by 
—X, leaves the interaction term —F(t)x in the Hamilto- 
nian ^ invariant, so that a phase shift by tt is equivalent 
to spatial inversion. Needless to say, one can also select 
other values of if and produce results not predictable by 
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FIG. 10: (Color online) As Fig. H but for ip = Tv. 



simple symmetry considerations. For example, the out- 
come of a passage through an avoided crossing may de- 
pend on the particular instantaneous phase at which this 
anticrossing is met. More generally, the possible effects 
of the phase of the carrier frequency with respect to the 
pulse envelope closely resemble corresponding effects en- 
countered in laser-atom interaction 1461. 



V. CONCLUSIONS AND FURTHER VISIONS 

In this paper we have discussed mechanisms which gov- 
ern the response of a single particle in a cosine lattice to 
pulsed homogeneous forcing, with a view toward con- 
trolling weakly interacting Bose-Einstein condensates in 
shaken optical lattices. While several successful experi- 
ments with ultracold atoms in strongly ac-forced optical 
lattices have already been reported 041, IHI, [13, EO, Ell , 
the systematic exploration of the possibilities of coherent 
control opened up, e.g., by deliberate pulse shaping, is 
likely to break new ground in ultracold-atom physics. 

This optimistic view is suggested by a simple parallel: 
Atomic and molecular physics in itself, not invoking the 
use of lasers, already is a fascinating and important sub- 
ject, but it becomes infinitely more rich when lasers come 
into play, allowing one, on the one hand, to do precision 
spectroscopy, and to perform deliberate state manipula- 
tions on the other. By the same token, ultracold atoms in 
optical lattices offer access to much fundamental physics, 
but there is a host of further options when applying in- 
ertial forces, either in the form of an ac drive for creating 
dressed matter waves, or in the form of carefully designed 
pulses in order to exert active coherent control. 
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Here we have focused on elementary mechanisms of 
control deriving from pulses with a smooth envelope, 
leading to adiabatic motion of a wave packet's momen- 
tum distribution on quasienergy surfaces created by spa- 
tioteniporal Bloch waves and to deviations from adia- 
baticity which can be purposefully exploited for reaching 
target states which may not be accessible by other means. 
The accompanying Floquet picture offers the distinct ad- 
vantage that it enables one to adapt many concepts devel- 
oped for the theoretical description of laser-matter inter- 
action, such as the notion of tt pulses. In our opinion, the 
actual implementation and observation of such tt pulses 
with dilute Bose-Einstein condensates in shaken optical 
lattices constitutes an experimentum crucis: If this can 
be done, many further related control scenarios will be 
equally viable, also involving variations of the driving 
frequency. 

As a future perspective it seems particularly rewarding 
to also carry over advanced strategies devised for control- 
ling molecular dynamics and even chemical reactions by 
specifically engineered laser pulses [25l - [29| . Such tech- 
niques often involve feedback loops for optimizing the 
pulse characteristics with the help of genetic learning al- 
gorithms; this approach is tantamount to "teaching lasers 



to control molecules" [1^. When working with Bose- 
Einstein condensates in optical lattices, control pulses 
can be applied by piezo-electrically juggling the posi- 
tion of a mirror which refiects the lattice-generating laser 
beam back into itself, thus giving rise to a moving stand- 
ing light wave in the laboratory frame, and to a corre- 
sponding inertial force in the comoving frame of refer- 
ence. By sheer analogy, one then would "teach mirrors 
to control condensates." In this manner one could assess 
the "reachability problem" for condensates in optical lat- 
tices and explore whether one can reach any preselected 
final-state distribution with the help of a proper pulse se- 
quence. In particular, it would be of major interest to ap- 
ply such strategies for creating mesoscopic Schrodinger- 
cat-like states, that is, quantum superposition states of 
condensates. 
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